This is where the final project report write-up goes.

Before you submit, make sure everything runs as expected.

You can add sections as you see fit. Make sure you have a section called “Introduction” at the beginning and a section called “Conclusion” at the end. The rest is up to you!

##Introduction - Load the tidyverse, ggplot, and rtweet packages

library(tidyverse)
library(ggplot2)
library(rtweet)
library(readr)

This data set was scraped from WineEnthusiast, a website that reviews and rates many differet types of wines.

wines <- read.csv(file = '../data/winemag-data-130k-v2.csv')[,-1]
set.seed(19630217)
wine_sample<- sample_n(wines, 1000)

EDA (correlation priceXpoints, with DataExplorer library? using (this)[https://datascienceplus.com/blazing-fast-eda-in-r-with-dataexplorer/])

wines %>% 
  ggplot() +
    geom_point(mapping = (aes(x = points, y = price)), na.rm = T)

wines %>%
      summarize(mean(price, na.rm=TRUE), 
                min(price, na.rm=TRUE),
                max(price,na.rm=TRUE), 
                sd(price, na.rm=TRUE))
wines %>%
      summarize(mean(points, na.rm=TRUE), 
                min(points, na.rm=TRUE),
                max(points,na.rm=TRUE), 
                sd(points, na.rm=TRUE))

Select the provinces based on points and Select the best province for wine based on the average points of the sample size.

#find the average number of points across the 1,000 samples

wine_per_province <- wine_sample %>% 
  select(province, points) %>% 
  summarise(points = mean(points))
wine_per_province

#Find the best province for wine using the average points across the 1,000 samples

best_province <- wine_sample %>% 
  group_by(province, points) %>% 
  filter(points > 88.669)
best_province  

Rating distribution

Bang for your buck, by variety

#sort by price, then points
#want to do an interaction variableor somethin?
#wine_cheap_but_good <- 
wines %>% 
  group_by(variety) %>% 
  summarise(mean(points))# %>% 
#  arrange(points)
  
  

##Conclusion

LS0tDQp0aXRsZTogIkZpbmFsIFJlcG9ydCBmb3IgSW50cm8gdG8gRGF0YSBTY2llbmNlIg0Kb3V0cHV0OiBodG1sX25vdGVib29rDQotLS0NClRoaXMgaXMgd2hlcmUgdGhlIGZpbmFsIHByb2plY3QgcmVwb3J0IHdyaXRlLXVwIGdvZXMuIA0KDQpCZWZvcmUgeW91IHN1Ym1pdCwgbWFrZSBzdXJlIGV2ZXJ5dGhpbmcgcnVucyBhcyBleHBlY3RlZC4NCg0KWW91IGNhbiBhZGQgc2VjdGlvbnMgYXMgeW91IHNlZSBmaXQuIE1ha2Ugc3VyZSB5b3UgaGF2ZSBhIHNlY3Rpb24gY2FsbGVkICJJbnRyb2R1Y3Rpb24iIGF0IHRoZSBiZWdpbm5pbmcgYW5kIGEgc2VjdGlvbiBjYWxsZWQgIkNvbmNsdXNpb24iIGF0IHRoZSBlbmQuIFRoZSByZXN0IGlzIHVwIHRvIHlvdSENCg0KDQoNCiMjSW50cm9kdWN0aW9uDQotIExvYWQgdGhlIGB0aWR5dmVyc2UsIGdncGxvdCwgYW5kIHJ0d2VldGAgcGFja2FnZXMNCmBgYHtyLCBtZXNzYWdlPUZBTFNFLCB3YXJuaW5nPUZBTFNFfQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KGdncGxvdDIpDQpsaWJyYXJ5KHJ0d2VldCkNCmxpYnJhcnkocmVhZHIpDQpgYGANCg0KDQpUaGlzIGRhdGEgc2V0IHdhcyBzY3JhcGVkIGZyb20gV2luZUVudGh1c2lhc3QsIGEgd2Vic2l0ZSB0aGF0IHJldmlld3MgYW5kIHJhdGVzIG1hbnkgZGlmZmVyZXQgdHlwZXMgb2Ygd2luZXMuICANCg0KLSBUaGlzIGRhdGFzZXQgaW5jbHVkZXMgaW5mb3JtYXRpb24gb2Ygb2YgMTMwLDAwMCB3aW5lIHJldmlld3Mgd2l0aCAxMCBkaWZmZXJlbnQgZGF0YSBmaWVsZHMuIA0KICANCg0KYGBge3J9DQp3aW5lcyA8LSByZWFkLmNzdihmaWxlID0gJy4uL2RhdGEvd2luZW1hZy1kYXRhLTEzMGstdjIuY3N2JylbLC0xXQ0KYGBgDQoNCmBgYHtyfQ0Kc2V0LnNlZWQoMTk2MzAyMTcpDQp3aW5lX3NhbXBsZTwtIHNhbXBsZV9uKHdpbmVzLCAxMDAwKQ0KYGBgDQoNCkVEQSAoY29ycmVsYXRpb24gcHJpY2VYcG9pbnRzLCB3aXRoIGBgYERhdGFFeHBsb3JlcmBgYCBsaWJyYXJ5PyB1c2luZyAodGhpcylbaHR0cHM6Ly9kYXRhc2NpZW5jZXBsdXMuY29tL2JsYXppbmctZmFzdC1lZGEtaW4tci13aXRoLWRhdGFleHBsb3Jlci9dKQ0KYGBge3J9DQp3aW5lcyAlPiUgDQogIGdncGxvdCgpICsNCiAgICBnZW9tX3BvaW50KG1hcHBpbmcgPSAoYWVzKHggPSBwb2ludHMsIHkgPSBwcmljZSkpLCBuYS5ybSA9IFQpDQpgYGANCg0KYGBge3J9DQp3aW5lcyAlPiUNCiAgICBzdW1tYXJpemUobWVhbihwcmljZSwgbmEucm09VFJVRSksIA0KICAgICAgICAgICAgICBtaW4ocHJpY2UsIG5hLnJtPVRSVUUpLA0KICAgICAgICAgICAgICBtYXgocHJpY2UsbmEucm09VFJVRSksIA0KICAgICAgICAgICAgICBzZChwcmljZSwgbmEucm09VFJVRSkpDQpgYGANCmBgYHtyfQ0Kd2luZXMgJT4lDQogICAgc3VtbWFyaXplKG1lYW4ocG9pbnRzLCBuYS5ybT1UUlVFKSwgDQogICAgICAgICAgICAgIG1pbihwb2ludHMsIG5hLnJtPVRSVUUpLA0KICAgICAgICAgICAgICBtYXgocG9pbnRzLG5hLnJtPVRSVUUpLCANCiAgICAgICAgICAgICAgc2QocG9pbnRzLCBuYS5ybT1UUlVFKSkNCmBgYA0KDQpTZWxlY3QgdGhlIHByb3ZpbmNlcyBiYXNlZCBvbiBwb2ludHMgIGFuZCBTZWxlY3QgdGhlIGJlc3QgcHJvdmluY2UgZm9yIHdpbmUgYmFzZWQgb24gdGhlIGF2ZXJhZ2UgcG9pbnRzIG9mIHRoZSBzYW1wbGUgc2l6ZS4gDQoNCiNmaW5kIHRoZSBhdmVyYWdlIG51bWJlciBvZiBwb2ludHMgYWNyb3NzIHRoZSAxLDAwMCBzYW1wbGVzDQpgYGB7cn0NCndpbmVfcGVyX3Byb3ZpbmNlIDwtIHdpbmVfc2FtcGxlICU+JSANCiAgc2VsZWN0KHByb3ZpbmNlLCBwb2ludHMpICU+JSANCiAgc3VtbWFyaXNlKHBvaW50cyA9IG1lYW4ocG9pbnRzKSkNCndpbmVfcGVyX3Byb3ZpbmNlDQpgYGANCg0KDQojRmluZCB0aGUgYmVzdCBwcm92aW5jZSBmb3Igd2luZSB1c2luZyB0aGUgYXZlcmFnZSBwb2ludHMgYWNyb3NzIHRoZSAxLDAwMCBzYW1wbGVzDQpgYGB7cn0NCmJlc3RfcHJvdmluY2UgPC0gd2luZV9zYW1wbGUgJT4lIA0KICBncm91cF9ieShwcm92aW5jZSwgcG9pbnRzKSAlPiUgDQogIGZpbHRlcihwb2ludHMgPiA4OC42NjkpDQpiZXN0X3Byb3ZpbmNlICANCmBgYA0KDQoNClJhdGluZyBkaXN0cmlidXRpb24NCg0KYGBge3J9DQoNCmBgYA0KDQpCYW5nIGZvciB5b3VyIGJ1Y2ssIGJ5IHZhcmlldHkNCmBgYHtyfQ0KI3NvcnQgYnkgcHJpY2UsIHRoZW4gcG9pbnRzDQojd2FudCB0byBkbyBhbiBpbnRlcmFjdGlvbiB2YXJpYWJsZW9yIHNvbWV0aGluPw0KI3dpbmVfY2hlYXBfYnV0X2dvb2QgPC0gDQp3aW5lcyAlPiUgDQogIGdyb3VwX2J5KHZhcmlldHkpICU+JSANCiAgc3VtbWFyaXNlKG1lYW4ocG9pbnRzKSkjICU+JSANCiMgIGFycmFuZ2UocG9pbnRzKQ0KICANCiAgDQpgYGANCg0KDQoNCg0KIyNDb25jbHVzaW9uDQo=